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Abstract. We study thermal convection in a colloidal glass of Laponite in formation. Low concentration 
preparation are submitted to destabilizing vertical temperature gradient, and present a gradual transition 
from a turbulent convective state to a steady conductive state as their viscosity increases. The time spent 
under convection is found to depend strongly on sample concentration, decreasing exponentially with mass 
fraction of colloidal particles. Moreover, at fixed concentration, it also depends slightly on the pattern 
selected by the Rayleigh Benard instability: more rolls maintain the convection state longer. This behavior 
can be interpreted with recent theoretical approaches of soft glassy material rheology. 

PACS. 44.25.+f Natural convection. - 61.43. Fs Glasses. - 61.20. Lc Time-dependent properties; relaxation. 



1 Introduction 



Many systems in nature, such as glasses, spin-glasses, col- 
loids and granular materials, present a very slow relax- 
ation towards equilibrium and any physical property of 
these materials evolves as a function of time: they are ag- 
ing. For example, when a glassy system is quenched from 
above its glass transition temperature T g to a tempera- 
ture lower than T g , any response function of the material 
depends on the time t w elapsed from the quench pQ. An- 
other example of aging is given by colloidal glasses, whose 
properties evolve during the sol-gel transition which may 
last several months [2]. 

Although great experimental and theoretical efforts 
have been paid for years to understand the physics of 
glassy systems, a global coherent picture is still laking. 
Some recent theories [3] based on the description of spin 
glasses by a mean field approach proposed to extend the 
concept of thermodynamic temperature using a Fluctu- 
ation Dissipation Relation (FDR) which generalizes the 
Fluctuation Dissipation Theorem (FDT) for a weakly out 
of equilibrium system (for a review see ref. |4I5I6| V This 
novel approach has been quite successful, and is already il- 
lustrated by a few experiments in real materials |7l8l9ll0fTT| 
where FDR could be measured and compared to theoret- 
ical predictions. The main difficulty of these experiments 
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is the study of intrinsic fluctuations, since no long-time 
averaging can be done in an aging system. 

It has been proposed that when a gentle forcing is ap- 
plied to these systems, they could reach a stationary state 
where no aging can be observed, although they are still 
out of equilibrium [3112] . In this state, the fluctuation- 
dissipation ratio would be similar to that of the glassy 
phase, with a two temperatures scenario in the case of 
p-spins glasses for example |12j. This theoretical idea is 
very appealing to the experimentalist as it would ease the 
fluctuation measurements, allowing long time averaging 
in a stationary system. The drawback of typical methods 
is that forcing the system, by adding a shear force for 
example, induces great perturbations on the fluctuation 
measurements due to unwanted vibrations, electrical dis- 
turbances, etc. inevitably produced by the forcing method. 

Thermally driven flows may present a bypass for these 
problems: the classic Rayleigh Benard instability for ex- 
ample, can be used as the forcing method. Heating the 
sample from below, one triggers convective motion, and 
therefore a small shear associated to the flow [13]. The 
advantage of this method is that it does not imply any 
moving parts or motors, which would inevitably produce 
additional fluctuations. 

In the study of soft glassy material, Laponite [H] has 
often been use as a standard example of colloidal glass, 
presenting many similarities with structural glass aging 
0. Rheology |9ll5ll6ll7j . light and X-Ray diffusion |2ll5ll6li8] . 
electrical measurements [819] have been used to character- 
ize aging of this material, and evidenced many interesting 
aspect of its glassy behavior. Nevertheless, no attention 
has been paid up to now to its thermal properties or to 
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its behavior in a convective flow. We propose here a in- 
vestigation in this direction, with the following questions 
in mind: are thermal properties of Laponite aging too ? 
Can convective flow modify (delay, accelerate, stop) aging 
behavior ? Would it make a suitable driving method to 
study the fluctuation-dissipation ratio of a forced glassy 
system ? 

This article is organized as follow: we first describe 
the experimental set-up (Sect. [2]), then study the basic 
properties relevant for convection in Laponite (Sect. [3|). 
Sect. U details the experimental results of convection in 
Laponite, before the concluding section (Sect. [5]). 



instruments are remote controlled with a computer via 
GPIB interface, which allows sampling temperature of the 
various sensors at 0.1 Hz and controlling hot temperature 
with a PID algorithm: we achieve better than 0.5 K tem- 
perature uniformity and 0.03 K stability. The temperature 
of the external environment is stabilized to better than 
0.5 K over 24 h. To minimize heat exchange of the cell 
with the environment, the hot and cold plate temperature 
are set symmetrically to the room one, and both plates 
are imbedded in a tick polystyrene insulation. In the worst 
case (conductive configuration) , heat losses represent 30 % 
of total heating power. 



2 Experimental setup 



2.2 Pattern visualization 



Cold plate (water circuialion 1 



Temperature sensors (PtlOO) 




Hot plate (heater foil) 



Fig. 1. Experimental setup: the convection cell is rectangular 
with aspect ratios close to 6 : 1 and 3 : 1 (L = 145 mm, 
w — 68 mm, h = 23.5 mm). The cold plate temperature is 
regulated with a water circulation, while the other plate is 
heated with a flat sheet resistance. Several PtlOO sensors are 
used to measure temperature stability and uniformity. 



2.1 Convection cell 

The experimental setup is illustrated in FigJTJ The con- 
vection cell is rectangular: (Lxwxh) 145 mm long, 68 mm 
deep and 23.5 mm high. The lateral walls are made of 
10 mm thick plexiglas, and top/bottom plates are made 
of 16 mm thick stainless steeo Several platinum sensors 
(RTD-100) are placed in the two metal plate bodies, close 
to the cell surface, to measure temperature gradient and 
uniformity. All measurements are made with a Keithlcy 
2000 multimeter and scan card in 4 wires configuration. 
The cold plate temperature is controlled with a circula- 
tion from a water bath (Nestlab RTE-210), allowing bet- 
ter than 0.1K temperature uniformity and stability. The 
other plate is heated with a 18 W flat sheet resistance of 
matched size. This resistance is connected to a DC power 
supply (HP E3632A) in a 4 wires configuration, which 
allows precise heating power setting and monitoring. All 



1 Although stainless steel has a poor thermal conductivity 
with respect to the other metals usually used in a convection 
cell, it is still more than one order of magnitude more conduc- 
tive than water and ensures chemical neutrality in Laponite 
preparation aging. 



To perform a visualization of the convective flow pattern, 
we propose here a very simple technique. We produce a 
horizontal laser light sheet using a cylindrical lens and a 
diode pumped solid state laser (Mcllcs Griot 58GCS405, 
5mW output at 532 nm). This light sheet crosses the cell 
along its shorter horizontal length, close to the bottom 
plate (at z = 2 mm). The temperature gradient inside the 
cell induces an optical index gradient, mainly in the (x,z) 
plane, deviating the light from its straight path. The light 
sheet is thus modulated by the temperature pattern, and 
its projection on a screen can be monitored with a camera: 
the line obtained represent the average of the temperature 
gradient in the y direction (shorter horizontal length of the 
cell). The slight horizontal deviation cannot be measured 
since it takes place inside the light plane, but the verti- 
cal deviation is easily recorded: for a typical temperature 
difference AT = 8K in the cell filled with water, com- 
puted and observed deviation is 2 x 10 -3 rad, when beam 
pointing stability is better than 2 x 10~ 5 rad. Thus rela- 
tive variation up to 1 % can be detected on the vertical 
temperature gradient d z T. 

This simple technique is sufficient to describe the flow 
pattern main characteristics such as size and periodicity of 
spatial structures, time evolution, etc. Quantitative mea- 
surement of the vertical temperature gradient can even be 
made when the deviation is compared to reference situa- 
tions (zero and uniform temperature gradient). Eventu- 
ally, performing a systematic vertical scanning would al- 
low the complete determination of temperature field, but 
this technique is beyond the scope of the present work. 
We present in figure [5] an example of this visualization 
technique in the case of 6 steady convection rolls. 



2.3 Laponite preparation 

A great attention is paid to prepare the Laponite sam- 
ples always in the same way, minimizing the source of 
non reproducibility. Preparation are done in a clean N2 
atmosphere, to avoid CO2 contamination. Laponite RD 
powder is mixed in a given mass fraction with distillate 
water for 1 h under strong magnetic stirring. Concentra- 
tion C ranges from (1.8 to 2.7) wt% in these experiments, 
leading to the slow formation of a colloidal glass in times 
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Fig. 2. Pattern visualization: A horizontal laser light sheet crosses the cell along its shorter horizontal length, 2 mm above the 
bottom plate. It is deviated by the temperature gradient integrated along its path through the cell. Its projection on a slanted 
screen (to magnify deviation) follow the temperature pattern, and can be recorded with a digital camera. We extract from the 
raw image the position of the maximum intensity line (thin red line in the center of the light sheet image). After appropriate 
rescaling, we can plot on the right figure the vertical temperature gradient (normalized by the measurement corresponding to 
the conductive case) along the cell. We present here the example of 6 regular convection rolls. 



ranging from months to hours. The solution, initially fairly 
liquid, is injected in the cell through a 1 /jm pore filter 
(Whatman Puradisc 25 GD, 1 /xm GMF-150) under 2 bar 
pressure. This event sets the initial time of the experiment. 
Filling and setting up operations take about 5 min, which 
is considered as the time accuracy of our measurement. 
To minimize thermalization time of the hot (cold) plate, 
it is preheated (precooled) to the desired temperature be- 
fore filling. Thanks to the symmetry of temperature limits 
around the room one, the mean temperature of the solu- 
tion does not change and its volume is constant. 



3 Physical properties of Laponite samples 

In a first set of experiments, we studied the time evo- 
lution of some physical properties relevant in Rayleigh 
Bcnard convection of Laponite preparations, in order to 
see if any aging behavior can be evidenced in these ob- 
servables of the system, as for its viscous or electrical 
properties. First of all, we measured the thermal con- 
ductivity A and the thermal diffusivity k using our ex- 
perimental setup in a conduction configuration: the cell 
is flipped upside down so that heating takes place from 
above, resulting in a vertical stabilizing temperature gra- 
dient, thus a conductive regime. The heat flux is imposed 
in a series of steps and the temperature difference across 
the cell is recorded for a 2.2 wt% preparation of Lapo- 
nite and for pure water as a calibration measurement. The 
transients in AT let us estimate the thermal diffusivity, 
while the stationary states lead to the conductivity. But 
as can be expected from their composition, none of the 
samples can be distinguished from pure water: no aging 
is observed and we measure A = (0.6 ± 0.02) W-m^-K" 1 
and K = (1.41 ±0.07) x 10~ 7 m 2 -s" 1 at 25 °C. 

Similarly, the thermal expansion coefficient a and den- 
sity p of Laponite preparation don't age and are equal to 
that of water within a few percent: for a 2.2 wt% prepa- 



ration at 25 °C, we measured p = (1 010 ± 3) kg-m 3 and 
a = (2.0 ±0.1) x lO^K" 1 . 

For a ideal newtonian fluid of kincmanic viscosity 
the control parameter of the Rayleigh Bcnard instability 
is the Rayleigh number, defined as 



Ra = 



gaATh 3 



(1) 



where g is the gravity constant and h the height of the con- 
vection cell [13]. For our samples, the rheological behavior 
would thus be the only parameter presenting substantial 
aging and be responsible for the convective properties of 
Laponite preparations: convection should be present as 
long as the Rayleigh number is greater than the critical 
threshold Ra c , that is as long as the viscosity is lower than 
the corresponding critical viscosity v c . Using Laponite pre- 
paration properties, AT — 8 K and assuming Ra c ss 2000 
[To] , we compute v c ps 7.5 x 10 _4 m 2 -s~ 1 . 

We present in figure [3] a measurement of the complex 
viscosity v* of a 2.2 wt% preparation at 25 °C. These re- 
sults are obtained with a commercial Bohlin rheometer 
(C-VOR 150) in cone plane configuration (60 mm diam- 
eter for 2° aperture) , under a imposed strain of 5 % and 
different frequencies / ranging from 0.1 Hz to 10 Hz. In 
this configuration, we test the linear viscoelastic response 
of the sample to harmonic shearing, and measure its elas- 
tic (G") and viscous (G") moduli. The complex viscosity 
v* is then defined as 



1 VG" 2 + G" 2 

P 2^7 



(2) 



For an ideal Newtonian liquid, this complex viscosity is in- 
dependent of frequency and equal to the classic viscosity, 
whereas for a Hook solid, it scales as the inverse of fre- 
quency The aging of Laponite rheological properties can 
be seen as a transition from a liquid like state (the vis- 
cosity being roughly independent of frequency and gradu- 
ally increasing) to a solid like state, where v* decreases as 
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Fig. 3. Complex viscosity of a 2.2 wt% Laponite preparation, v* is plotted as a function of time for various frequencies (left), 
and as a function of frequency for different aging times (right). We observe a transition between an initial liquid state, where 
v* is roughly independent of /, and a solid like behavior for larger times, where v* scale as J -0 - 85 (dashed line). Convection 
should occur as long as v* is smaller than v c , thus about 15 h to 16 h at this specific concentration. Low time / high frequency 
data is not available due to inertia limitations of the measurement. 



[XT] As seen on figure [31 this transition occurs for 
aging times around 15 h to 16 h for a 2.2 wt% preparation, 
when v* abruptly crosses the v c threshold at lower fre- 
quencies. As the apparition of the viscoclastic properties 
of Laponite preparations occurs lately, when the viscosity 
is already quite high, the convective behavior of the sam- 
ple can be considered as that of a casual liquid whose vis- 
cosity increases with times. The elastic component might 
only play a role in the last part of the experiment, just 
before the onset of the conductive state. 



4 Convection with Laponite 

To study the influence on convection of the aging of La- 
ponite preparations, let us focus on the time evolution of 
the heat transfer in a convective configuration: we heat the 
cell from below and impose a constant temperature differ- 
ence AT = 8 K, and record the injected power. A simple 
ratio with the equivalent quantity in a conduction con- 
figuration (with necessary corrections to compensate for 
heat losses) define the Nusselt number Nu of this sample. 
In figure [U we present its time evolution for a few La- 
ponite preparations with concentrations between 2.0 wt% 
and 2.4 wt%. Initially, the solution has a very low vis- 
cosity (sligthly higher than water) and the heat transfer 
is efficient, leading to a Nusselt number close to the one 
reached with water in equivalent conditions: Nu = 7.4 
(for Ra = 1.5 x 10 6 ). As the sample ages, its viscosity in- 
creases, reducing the Raylcigh number and thus the Nus- 
selt number, up to the moment where convection disap- 
pear: the Rayleigh number is smaller than the critical one, 
and the heat transfer is purely conductive (Nu — 1 ). This 
instant define a characteristic time: the duration of convec- 
tion for a specific sample will be denoted as t c . Although 
the concave shapes of all the curves are quite similar, no 




30 40 
Time, t/h 

Fig. 4. Aging of the Nusselt number: samples with different 
concentration C (in wt%) are submitted to a AT — 8K tem- 
perature difference in convective configuration, and we plot 
the time evolution of the Nusselt number. Initially the heat 
transfer is efficient, close to the one reached with water in sim- 
ilar conditions, but the increasing viscosity of Laponite slows 
down the convective motion to eventually stop the Rayleigh 
Benard instability for a given time t c (C). In the inset, time 
axis is rescaled by t c for each ploted concentration (from top 
to bottom: C = (1.8, 2.05, 2.2 and2.6)wt%). Although the 
curves are very similar, we notice a slight decreasing concavity 
of Nu(t/t c ) with C. 



simple rescaling was found in order to get a master curve: 
in the inset of figure [4j we present a simple rcnormaliza- 
tion of the time axis by t c , evidencing the slight decreasing 
concavity of Nu(t/t c ) with the concentration C. 

In figure Owe plot the convection time t c as a function 
of sample concentration C, for AT = 8K. t c turns out 
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Fig. 5. Convection time t c as a function of sample concen- 
tration C. Measurements are all done with AT = 8K.The 
vertical logarithmic scale underline that t c is a simple expo- 



nential function of C: 
Co = 0.17 wt%. 



t c = t e~ c/Co with to = 7.8 x 10 6 h and 



to be a simple exponential function of C: a best fitting 
procedure leads to 



t e 



-ci c 



(3) 



with t = 7.8 x 10 6 h and C„ = 0.17 wt%. This behavior of 
typical aging time versus concentration is consistent with 
previous observations on Laponite aging [2], although co- 
efficients t and Co are quite different since the physical 
properties studied are very different. In order to perform 
several experiments in a reasonable time, but with a typ- 
ical aging time much larger than the main thermal dif- 
fusion scale t = h 2 /n, we choose to work with 2.2 wt% 
samples, for which t c w 20 h. This way, we can make one 
experiment per day and still have more than one order of 
magnitude between t c and r = 1.1 h. 

If we repeat this same simple convection experiment 
with several samples of concentration C = 2.2 wt%, we 
see an important dispersion in t c : as shown in figure [6l t c 
varies between 18 h and 22 h. As emphasized in the para- 
graph on Laponite preparation, we take great care to make 
every preparation in the same way, and the differences be- 
tween these sample history arc tiny. The precision on con- 
centration for example, better than 0.002 wt%, exclude a 
interpretation of this dispersion in terms of sample compo- 
sition : indeed, this imprecision leads to a 0.5 h incertitude 
according to eq[3j Moreover, we see that all the curves ex- 
actly superpose in a first regime, for the first 5 h, and only 
then the behavior of each experiment changes to fall more 
or less quickly into the conductive state. The visualization 
of the pattern gives an interesting insight of the mecha- 
nism at work: the convective flow is turbulent initially, 
with small structures quickly evolving in time (the typical 
time of decorrelation of the pattern is 5 s after a 2 h aging). 
This state is reproducible among samples and leads to an 
efficient heat transfer. As the viscosity of the sample in- 




Fig. 6. Reproducibility of experiments: time evolution of Nus- 
selt number, for AT = 8 K and C = 2.2 wt%. As shown by the 
visualization technique, the convective flow is turbulent when 
all curves initially superpose. The final decaying of Nu to 1 is 
pattern dependent: t c is smaller for fewer rolls. 



creases, the flow slowly evolves towards a more structured 
pattern, involving steady like convection rolls. The pre- 
ferred pattern in a cell with 6 : 1 aspect ratio should be 
of 6 steady convection rolls, but the lateral equilibration 
time of the cell is about 42 h, longer than the existence 
time of the pattern! Therefore, other patterns may be ini- 
tially selected by the turbulent flow, and can survive till 
the end of convection. This is what is observed in the last 
hours of convection of our samples: some patterns with 
5,6 or 7 rolls have been observed, as well as other more 
complicated structures (generally unsteady - not shown). 
It turns out qualitatively that patterns with less rolls, fall 
quicker in a conductive state. The interpretation of the 
dispersion of the convection times t c would be the follow- 
ing: when smaller flow structures are implied, the average 
shear rate of the flow is higher, thus delaying the normal 
aging of the sample and maintaining convection for longer 
times. 



5 Discussions and conclusions 

We present a study of low concentration Laponite prepa- 
rations submitted to a destabilizing vertical temperature 
gradient. Except for rheological behavior, all the relevant 
properties of the samples are those of water, which in 
similar experimental conditions present turbulent convec- 
tion (Ra = 1.5 x 10 6 , Nu = 7.5). As measured with a 
standard rheometer, the complex viscosity of the samples 
evolves in time from fluid like to a solid like behavior. The 
onset of viscoelastic properties coincide with the increase 
of viscosity beyond the critical value for which convection 
disappear. 
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In the convection cell, this transition can be followed 
with the decaying of the Nusselt number from 6 to 1, as the 
heat transfer relaxes from a turbulent convective flow to 
a purely conductive state. The duration of the convective 
state t Cl depends strongly on the sample concentration : it 
decreases exponentially on more than 2 orders of magni- 
tude when C changes between 1.8 wt% and 2.7wt%. This 
sensitivity shows that convection is a very sensitive tool 
to investigate Laponite aging. 



0.5 




25 



Time, t/h 

Fig. 7. The mean shear rate to which the sample is submitted 
to is estimated from eq.|4](see text for details). Although quite 
small, this small forcing associated to the thermally driven 
flow delays the aging of the sample, that is the formation of 
the colloidal glass. 

For a given concentration, t c is longer than expected 
from rheological measurements : while v* crosses the crit- 
ical viscosity v c around 15 h, convection last about 20 h. 
Moreover, t c presents an important dispersion, which is 
attributed to the Rayleigh Benard instability pattern: con- 
vection can be up to 20 % longer when more rolls are im- 
plied. These effects can be interpreted considering the ef- 
fect of the flow on aging : the small shear associated to 
the thermally driven fluid motion can delay aging, de- 
laying the formation of the colloidal glass. Indeed, let us 
consider the global energy balance in the cell [13] : 

(|V«| 2 ) = ^ik(JV«-l) (4) 

The left term of this equation represent the space average 
of the velocity gradient, that is the mean shear rate to 
which the sample is submitted. Using the measurement of 
v* at 1 Hz (figure[3]) as an estimation for the kinematic vis- 
cosity, equation [1] to compute the Rayleigh number, and 
the average behavior of the Nusselt number of figure 
we can plot on figure [7] a rough estimation of the time 
evolution of the shear rate to which the sample is submit- 
ted. The order of magnitude of this shear rate is 0.3 s _1 . 
Although it is quite small, this forcing delays the aging 



of the preparation, and convection reveals this behavior 
by presenting a sensitivity to the flow pattern: more rolls 
(smaller spacial scale) present more shear for the sample, 
thus a longer convection state. 

In other experiments performed on Laponite prepara- 
tion to study their behavior under shear [5T], to get a 
stationary viscosity around 1 Pa-s, the shear rate needed 
is about 25 s -1 for concentrations of (3 to 3.5) wt%. The 
concentrations we probe in our experiment are much 
smaller, thus lowering the viscosity and probably the driv- 
ing needed to get a stationary state. Anyway, our shear 
rate stand 2 orders of magnitude below this value, and we 
couldn't reach such a stationary state. Moreover, increas- 
ing the convective shear to larger values appears to be 
quite difficult: according to eq. 5] it would mean increas- 
ing the Rayleigh and Nusselt numbers. Nu is known to 
scale roughly as Ra 1 ^ 3 cx ft. as long as Ra » Ra c , so that 
(|Vu| 2 ) is independent on the value of h, and the only way 
to increase it is to raise AT. We already use quite a high 
temperature difference (AT = 8K), and cannot pretend 
to gain an order of magnitude on this parameter with- 
out loosing the Boussinesq approximation |13j . Thus the 
only chance to reach a stationary convective state would 
be to lower even more the concentration, which implies 
according to cq.[3]that we would have to deal with month 
(year ?) like time scales... 

As a conclusion, let us eventually try to answer the 
few questions we had in mind when beginning this study. 
First, thermal properties of low concentration Laponite 
preparation are not aging, they are in fact those of water 
within a few percent. Second, convective flow can indeed 
modify the aging behavior, delaying the formation of the 
colloidal glass when more convection rolls are implied. But 
anyway, it probably won't make a suitable driving method 
to study the fluctuation-dissipation ratio of a forced glassy 
system: we are unable to reach a stationary state which 
would ease the measurement of fluctuations. 

However, this study gives a nice illustration of the 
coupling between flow and aging. It also sheds light on 
an interesting point: when forcing a soft glassy material, 
flow instabilities can play a significant role in the overall 
glassy behavior. In a wider field of view, this conclusion 
can be extended to many systems with inner processes 
(phase transitions, chemical reactions, etc.) occurring on 
time scales comparable to the flow ones. 
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